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We construct a quantum Ginsburg-Landau theory to study the quantum phases and transitions 
in electron hole bilayer system. We propose that in the dilute limit as distance is increased, there is 
a first order transition from the excitonic supcrfluid (ESF ) to the excitonic supersohd (ESS) driven 
by the collapsing of a roton minimum, then a 2nd order transition from the ESS to excitonic normal 
solid. We show the latter transition is in the same universality class of superfluid to Mott transition 
in a rigid lattice. We then study novel elementary low energy excitations inside the ESS. We find 
that there are two " supersolidon " longitudinal modes ( one upper branch and one lower branch ) 
inside the ESS, while the transverse mode in the ESS stays the same as that inside a ENS. We also 
work out various experimental signatures of these novel elementary excitations by evaluating the 
Debye- Waller factor, density- density correlation, specific heat and vortex -vertex interactions. For 
the meta-stable supersolid generated by photon pumping, we show that the angle resolved spectrum 
is dominated by the macroscopic super-radiance from its superfluid component, even it is just a very 
small percentage of the the whole system. This fact can be used to detect the metastable ESS state 
generated by photon pumping by a power spectrum experiment easily and without any ambiguity. 
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I. INTRODUCTION 

Recently, degenerate exciton systems have been pro- 
duced by different experimental groups with differ- 
ent methods in quasi-two-dimensional semiconductor 
GaAs/AlGaAs coupled quantum wells structurei'^i'^i'^i^. 
There are two different ways to generate the excitons. 
One is through photon pumping, then applying an elec- 
tric field along z direction to confine electrons in one 
quantum well and holes in another quantum wellii^i^. 
We call these kind of excitons as photon generated ex- 
citons. The second is through gate voltage^. We call 
these kind of excitons as gate voltage generated excitons. 
When the distance between the two quantum wells is suf- 
ficiently small, an electron in one well and a hole in the 
other well could pair to form an exciton which behaves 
as a boson at long distance. Now it was widely believed 
that the electron-hole bilayer (EHBL) is a very promising 
system to observe Bose-Einstein condensation (EEC) of 
excitons. For simplicity, we only discuss balanced el-hole 
bilayers where the electrons and holes have the same den- 
sity Ue = Uh, but may have different masses rUe ^ rrih- 
The case with different masses and with different den- 
sity is also very interesting and will be investigated in 
separate publications. 

There are two important dimcnsionless parameters in 
the EHBL. One is the dimcnsionless distance 7 = d/as ( 
as is the Bohr radius ) between the two layers. Another 
is rg which is the ratio of the kinetic energy over the po- 
tential energy in a single layer. The r^as is the typical 
interparticle distance in a single layer, ft is easy to see 
that the ratio of intralayer Coulomb Vn over the inter- 
layer Coulomb V12 interactions a = V'ii/Vi2 = d/rsCiB- 
So when the interlayer Coulomb interaction dominates 
a < 1, the EHBL is expected to exhibit the superfluid 
of excitons. If the density of excitons is sufhciently low ( 
large ), then the system is in a weakly coupled Wigner 



solid state at very large distance and become an BEC 
excitonic superfluid ( ESF ) at short distance ( Fig.l). 
An interesting problem is how the system evolves from 
the BEC ESF to the weakly coupled Wigner sohd as 
the distance increases. If an exciton is already formed, 
its kinetic energy K ~ m*(r ' potential energy 
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:(rjj,y-i - When K < P, namely, ^ < d/as, the 
EHBL could favor a excitonic ( or dipolar) normal solid 
(ENS) state. As argued above, when d/as < rg, the 
EHBL is in a ESF state. So in the intermediate distance 
< d/uB < Ts^, the system may favor a excitonic ( 
or dipolar ) supersolid (ESS) state. When d/as > rg, 
it will become the excitonic normal solid (ENS) due to 
the long range dipole-dipole 1/r^ repulsive interactions 
( Fig.l ). A hole system with ~ 20 — 30 was already 
realized^. If Ts ~ 20 3> 1 limit, there is a broad distance 
regime 4.5 < d/as < 20, the system could be in the ESS 
state. It becomes feasible to experimentally explore all 
the possible phases and phase transitions in the EHBL 
in the near future. 

It was originally pointed out in^^ that for any super- 
solid to exist, there must be quantum fluctuations gen- 
erated vacancies even in the ground state of a solid at 
T = 0. Obviously, the first candidate to search for su- 
persolid is near the phase boundary of superfluid '^He 
and solid '^He. The authors ini^ suggested that the su- 
persolid state leads to the non-classical rotational inertia 
(NCRI) observed in the torsional oscillator experiment, 
many other experiments such as neutron scattering, X- 
ray diffraction, mass flow, heat capacity, acoustic sound 
attenuation and so on are needed to uniquely distinguish 
the supersolid from other enormous number of much less 
interesting possibilities. A phenomenological quantum 
Ginzburg-Landau theoryi^ was developed to make pre- 
dictions on signatures of these experiments if the super- 
solid indeed exists in Helium 4. Unfortunately, so far. 
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all these experiments came as negative. It becomes in- 
teresting to see if a supersolid state can exist in other 
systems which also have both superfluid state and solid 
state. In this paper, we will point out that the excitons 
in EHBL maybe a very promising experimental system 
to observe the excitonic supersolid (ESS) near the phase 
boundary between the ESF and the ENS. Because the 
exciton's mass rriex = 0.37me is much smaller than that 
of '^He, so the zero point quantum fluctuations in the 
ENS in EHBL is even much larger than those in solid 
Helium 4. I will point out a new mechanism which is 
absent in '^He to generate repulsive excitonic vacancies 
in the ENS which leads to the intermediate excitonic su- 
persolid (ESS) phase in the dilute limit in the EHBL. 
Then we construct a quantum Ginzburg-Landau theory 
to study all the phases and quantum phase transitions 
in Fig.l. It is instructive to compare the EHBL system 
with the pseudo-spin sector in the bilayer quantum Hall 
system ( BLQH ) at total filling factor i^t — ^- Although 
the ESF in Fig.l shares some common properties with 
the corresponding ESF in the BLQlti^iiiii^, due to dif- 
ferent symmetries of the two systems, translational sym- 
metry breaking states are very different. For example, 
ESS and ENS phases which are the focus of this paper 
have no analogies in the BLQH. These crucial differences 
will be explicitly studied in the following. Furthermore, 
in sharp contrast to BLQH which is a stable system, for 
the photon generated excitons, the excitonic phases in 
Fig.l are just meta-stable states which will eventually 
decay by emitting lights. We will show that the angle 
resolved power spectrum from the internal photon can 
detect the ESS unanimously if it indeed exists. In fact, 
the characteristics of emitted photons is a very natural, 
feasible and unambiguous internal probe of all the three 
phases ESF, ESS and the ENS in Fig.l. 



The paper is organized as follows. In sec. II, we de- 
rive the quantum Ginsburg-Landau action to describe 
the transition from the ESF to the ENS driven by the 
collapsing of a roton minimum. Then In Sec. Ill, we will 
argue that in general there should be a ESS state in- 
tervening between the ESF and the ENS. In Sec. IV, by 
renormalization group analysis, we study the universal- 
ity class of zero temperature quantum phase transition 
from ENS to ESS driven by the distance. In Sec.V, we 
work out the elementary low energy excitations inside 
the supersolids. Then in the following sections, we study 
the experimental signatures of these low energy excita- 
tions by calculating the Debye- Waller factor in the X-ray 
scattering from the ESS in sec. VI, the density-density 
correlation function in the ESS in sec. VII. In Sec. VIII, 
we study the specific heat in the ESS. In Sec. IX, by 
performing a duality transformation to the vortex repre- 
sentation, we will study the vortices in the ESS. In Sec.X, 
we will present the photon emission pattern from the ESS 
formed by photon pumping generated excitons. Finally, 
we reach conclusions in Sec. XL 



II. THE ZERO TEMPERATURE TRANSITION 
FROM ESF TO ENS DRIVEN BY THE 
DISTANCE 

If ci(c2) is the electron annihilation operator in top ( 
bottom ) layer, then h2 = ^2 is the hole creation oper- 
ator in the bottom layer. The order parameter for the 
ESF is the p — h pairing ip{x,T) = (c|c2) = (cj/ij) = 
\/p + 5pe^^'^^''^\ The effective action inside the ESF is 
essentially the same as that in the pseudo-spin channel 
in BLQHiO; 

c[5p, e] = iSpdre + ipd(V0)2 + ]^5pVMSp (i) 

In the ESF state, it is convenient to integrate out 5p in 
favor of the phase field to get the phase representation: 

where the dispersion relation of the Goldstone modes in- 
cluding higher orders of momentum can be extracted: 

= [2pdVd{m^ (3) 

In the long wavelength limit, Vd((f ^ 0) ^ c ( c is 
a constant ) ( Fig.l ) leads to a capacitive term for the 
density fluctuation. The QMC calculationsli^ indeed find 
that there is a roton minimmn in the dispersion relation 
( Fig.l). 

Because the original instability comes from the density- 
density interaction Vd{q) , it is convenient to integrate 
out the phase field in favor of the density operator in the 
original action Eqn[TJ Neglecting the vortex excitations 
in 9 and integrating out the 9 in Eqn[l] leads to: 

^[M = \5p{-q, -^")[^ + VdimPiq.^n) (4) 

where we can identify the dynamic density-density cor- 
relation function: 

S{q,UJn) = {5p{-q,-UJn)5p{q^^n)) = 9 ^^'^2^/' \ 2 (^) 

+ V (q)q 

where v'^{q) ~ 2pdVd{q) is the spin wave velocity defined 
in EqnO 

From the pole of the dynamic density-density correla- 
tion function, we can identify the speed of sound wave 
which is exactly the same as the spin wave velocity. This 
should not be too surprising. As shown in liquid *He, 
the speed of sound is exactly the same as the phonon 
velocity. Here, in the context of excitonic superfluid, we 
explicitly prove that the sound speed is indeed the same 
as the spin wave velocity. 

From the analytical continuation iujn — > co + iS in 
EqnjSl we can identify the dynamic structure factor: 
S{q, uj) = S{q)S{uj — v{q)q) where S{q) — PdqT^/v{q) is the 
equal time density correlation function shown in Fig.l. 
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As g ^ 0, S{q) — > q. The Feymann relation in the ESF 
which relates the dispersion relation to the equal time 
structure factor is 



Q)(q) 



Lu{q) 



S{q) 



(6) 



which takes exactly the same form as the Feymann rela- 
tion in superfluid '^He. Obviously, the Vd{q) in the Fig. lb 
leads to the roton dispersion cj^ — q^Vd{q) in the Fig. la. 

Because the instability happens near q = qo instead of 
g = 0, so the transition in Fig.l is not driven by vortex 
unbinding transitions like in 3D XY model, so the vor- 
tices remain tightly bound near the transition. So inte- 
grating out the vortex excitations in 9 will only generate 
weak interactions among the density Sp: 



2"''^ ^' "''2pdq^ 
w{dpf + u{Sp)^ + ■ 



ydiq)]Sp{q,uj) 



(7) 



where the momentum and frequency conservation in the 
quartic and sixth order term is assumed. 

Expanding Vd{q) near the roton minimum qo leads to 
the quantum Ginsburg-Landau action to describe the 
ESF to the ENS transition: 
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C[Sp] - -Sp[ApUj'+r+ciq'^q'o) ]Sp-w{Spr+uiSp)^+- 



where Ar 



(8) 



„ .> which is non-critical across the tran- 



sition. In sharp contrast to the ESF to the pseudo-spin 
density wave (PSDW) transition in BLQH^", because of 
the lack of Z2 exchange symmetry between the two layers 
in EHBL, there is a cubic term in EqulS] It was explicitly 
shown ini^ that in the presence of both the cubic and the 
quartic terms, the favorable lattice is a triangular lattice 
instead of a square lattice in PSDW^*'. The generic tran- 
sition driven by the collapsing of roton minimum is from 
ESF to ENS instead of from the ESF to the excitonic 
supersolid ( ESS ). In the ESF, r)0, (ip) ^ 0, (6p) = 0, 
In the ENS, r < 0, {^j) = 0, (dp) = J2'g nge'^'^ where G 
is the 6 shortest reciprocal lattice vector of a triangular 
lattice ( Fig.l ). 



III. EXISTENCE OF ESS INTERVENING 
BETWEEN THE ESF AND ENS. 

In the section, I point out a mechanism to generate 
excitonic vacancies which lead to a narrow window of 
ESS intervening between the ESF and the ENS. As the 
distance increases to the critical distance dd, because 
the lattice constant VgaB is completely fixed by the pa- 
rameter which is independent of the distance which 
drives the transition, so the resulting state is likely to 
have vacancies with density ny{0) even at T = 0. By 
contrast, in solid Helium 4, the density is self-determined 



0)(q) 
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FIG. 1: The zero temperature phase diagram in the SEHB 
as the distance between the two layers increases. ESF where 
< ip >7^ 0, < >= stands for excitonic superfluid, ESS 
where < tp >^ 0, < >^ stands for exciton supersolid 
phase, ENS where < ip >— 0, < >^ stands for exci- 
ton normal solid phase, (a) Energy dispersion relation Lu{q) 
in these phases, (b) the bare dipole-dipole interaction Vd{q) 
in these phases. The cross in the ESS means the negative 
minimum value of Vd{q) is replaced by the ESS state. The 
order parameters are also shown. Also shown in far right in 
(a) is the equal time structure factor S{q). In fact, the insta- 
bility happens before the minimum touches zero. The generic 
transition driven by the collapsing of roton minimum is the 
ESF to ENS transition. However, as argued in the text, there 
could be a window of ESS intervening between the ESF and 
the ENS. 



by the pressure n 



dP 



TV, so the density and pres- 



sure go hand in hand, the solid "^He is likely to be com- 
mensurate. We expect that the vacancy-vacancy inter- 
action is also a repulsive dipole-dipole one. It is the 
condensation of these repulsively interacting vacancies at 
T ~ {) which leads to the SF mode ip{x, r) inside the in- 
commensurate ENS. This resulting state is the ESS state 
where {ip) ^ 0, {5p) ^ ( Fig.l ). As the distance in- 
creases to dc2 > dci,ny{0) — 0, the resulting state is a 
commensurate ENS whose lattice constant is still locked 
at rsQB (Fig.l). As distance increases further, the ENS 
will crossover to the two weakly coupled Wigner crystal. 
Note that because Vd((f ^ 0) — > c, the phonon spectrum 
is still w ~ q instead of ~ g^^^ as claimed in^. Very simi- 
lar argument was used ini^ to conclude that there should 
be zero-point quantum fluctuations generated vacancies 
in the pseudo-spin density wave ( PSDW ) in BLQH. By 
contrast, due to the absence of the cubic term, the lattice 
is a square lattice. In the PSDW, the vacancies are es- 
sentially fermionic holes, so can not condensate. It is the 
correlated hopping of vacancies in the active and passive 
layers in the PSDW state which leads to very large and 
temperature dependent drag consistent with the experi- 
mental data. 

The finite temperature phase diagram corresponding 
to Fig.l is shown in Fig. 2. At any finite T, the ESF will 
only have algebraic long range order and will turn into 
a normal liquid (NL) by a Kosterlitz-Thouless ( KT ) 
transition. ENS may go through a hexatic phase before 
melting into a NL phase. The finite temperature phases 
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FIG. 2: The finite temperature phase diagram in the SEHB 
as the distance between the two layers increases. ESF(A) 
means only algebraic off-diagonal long range order( ODLRO 
), ESS (A) means only algebraic ODLRO and algebraic trans- 
lational order. ENS(A) means only algebraic translational or- 
der. Hexatic phase has the algebraic orientational order, but 
no translational order. NL is the normal liquid phase. Tkt is 
the KT transition. Tv is the vector Coulomb gas transition. 
The dot is the zero temperature transition from the ENS to 
the ESS investigated in section 4. 



and phase transitions above the ESS phase in the inter- 
mediate distance dc\ < d < dc2 is also shown in Fig. 2. 



IV. THE ZERO TEMPERATURE TRANSITION 
FROM ESS TO ENS DRIVEN BY THE 
DISTANCE. 

The effective action to describe the ENS to ESS tran- 
sition at T = consistent with all the lattice symmetries 
and the global U{1) symmetry is: 



1 



(9) 



where r = d ~ dc2 ( Fig. 1), p„ is the normal den- 
sity, Uap — ^{daUp + dfjUa) is the linearized strain ten- 
sor, XaP'jS are the bare elastic constants dictated by 
the symmetry of the lattice, it has 2 independent elas- 
tic constants for a triangular lattice where XafSjS = 
XSafsS-ys + ^liSa-ySps + SasSpj) whcrc A and /i are Lame 
coefficients. In the ENS state r > 0, {ip{x,T)) = 0, the 
2 lattice phonon modes u(x, r) become the 2 ordinary 
ones. While inside the ESS state r < 0, {iIj{x,t)) ^ 0. If 
setting all the couplings between ip and vanish, the 
ijj sector describes the SF to Mott insulator transition 
in a rigid underlying two dimensional lattice studied 
inil. Under the Renormalization group ( RG ) trans- 
formation, r' — T/b^,x' = x/b and ip' = ip/Z. If we 
choose z — 2,Z = b~'^^'^, the g' = gb'^~''-. We also choose 
u'a — Ua/Z, then p'j = b~^pn, so the lattice phonon 
kinetic energy term is irrelevant near the QCP. It is 
easy to see g'q = b~'^/^~^ao, so Gq is always irrelevant. 
a'l = &^~''/^ai, so both g and ai's upper critical dimen- 
sion is (i„ = 2, so we can perform a e = 2 — d expansion 



in both g and ai. The RG equations are found to be: 
dg 2 '^'^1 ^ 



dl 



^g~cg ; 



dl 



(10) 



where c = 2mKdA'^~^ /h? is the same constant as that 
in the rigid model^^. So ai is exactly marginal with- 
out affecting the universality class. This is can also 
be understood as follows: because p„ is irrelevant, so 
we can simply integrate out Ua which only leads to a 
shift of the value of g. We conclude that the ENS to 
ESS transition stays in the same universality class of 
the superfluid to Mott insulator transition at c? = 2 
in a rigid lattice which has the mean field exponents 
z = 2,1^ = 1/2, = with logarithmic corrections. For 
example, the superfluid density inside the ESS should 
scale as ps \dc2 — d\^'^'^^~^'>'^ = \dc2 — d\ with a loga- 
rithmic correction. 

If neglecting the r dependence by setting Ua {x, r) = 
Ua{x),ipo{x,T) = '00 (a;), then EqnlH] reduces to the clas- 
sical action. For the classical case, x' = x/b^ip' = 
ii}/Z,u[^ = Ua/Z, if we choose Z = 6(2-'i)/2^ then 
g' = gb'^~'^,a'i — 6^~'^/^ai, so both g and ai's upper 
critical dimension is c?„ — 4. So in principle, a e = 4 — d 
expansion is possible for both g and ai , the putting e = 2 
for d — 2. It is known that due to the essential singu- 
larity of the KT transition, the specific heat exponent 
of the KT transition a — — oo < 0, the a^^ coupling 
is irrelevant, so the ENS to ESS transition remains to 
be Kosterlitz-Thouless (KT) transition at finite temper- 
ature. This conclusion is consistent with the RG analysis 
at r = in the last paragraph. The quantum critical 
scalings near the ENS to the ESS transition at d = dc2 
can be worked out along the similar lines in^S. 



V. THE LOW ENERGY EXCITATIONS IN THE 

ESS 

In this section, we will study the low energy elementary 
excitations in the ESS. Inside the ESS, {ijjo{x,T)) ~ a, 
we can write iPq{x,t) — -^/i + Spe^^^^''^'^ and plug it into 
the Eqnini Integrating out the massive magnitude 5p 
fluctuations and simplifying, we get the effective action 
describing the low energy modes inside the SS phase: 

+ \Hdref + pl,pda.9dp9]+a^pu^pidre (11) 

where k is the SF compressibility and is the SF 
stiffness which has the same symmetry as aP^p, ciap = 

'^ap + ^^^\p where S{^(k,LS) is the bare SF density cor- 
relation function. Obviously, the last term is the crucial 
coupling term which couples the lattice phonon modes 
to the SF mode. The factor of i is important in this 
coupling. By integration by parts, this term can also be 
written as aap{drUpdaO + drUadpO) which has the clear 
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physical meaning of the couphng between the SF velocity 
daO and the velocity of the lattice vibration d-rUp. It is 
this term which makes the low energy modes in the SS 
to have its own characteristics which could be detected 
by experiments. In this section, we neglect the topolog- 
ical vortex line excitations in Eqn illl In section IX, we 
will discuss these vortex line excitations in detail. In the 
following, we discuss triangular lattice specifically. 

For a triangular lattice, A^^^^ = \5ap5.^s + ^^{^af5ps + 
5as5pj) where A and u are Lame coefiicients, p% p — 

p'^^a,p,0'a,i3 = o-^a.p- In {Qi'^u) space, the Eqn[Tl] be- 
comes: 

^is = + (A + 2^)(7^]|u;(g,w„)P 

+ l[pn^l+ Pq^]\uM^n)\^ (12) 

where ui{q,uJn) = iqiUi{q,uJn)/q,Ut{q,uJn) = 



If setting a — 0, then Cj. reduces to the longitudinal 
phonon velocity vf^ = (A -f 2/i)/p„ and the superfluid 
velocity v'^ — Ps/ k respectively. Of course, the transverse 
phonon velocity v^p = p/pn is untouched. For notation 
simplicity, in the following, we just use Vp for vip. Inside 
the ESS, due to the very small superfluid density ps, it 
is expected that Vp > Vg. It is easy to show that v+ > 
Vp > Vs > f- and v'^ + v'^ > Vp + v'^, but v+v- ~ VpVs, 
so v+ +v- > Vp + Vg ( see Fig. 3 ). 

In principle, inelastic neutron scattering experiments 
or acoustic attenuation experiments can be used to detect 
the predicted the low energy excitation spectra in the SS 
shown in Fig. 3. 



VI. DEBYE- WALLER FACTOR IN THE X-RAY 
SCATTERING FROM THE ESS 

It is known that due to zero-point quantum motion in 
any NS at very low temperature, the X-ray scattering 
amplitude /(G) will be diminished by a Debye-Waller 
(DW) factor ~ e"^'^^^""^ at d = 2 where Ua is the lat- 
tice phonon modes in Eqn illl In Egn llH if the coupling 
between the u and 9 were absent, then the DW factor 
in the SS would be the same as that in the NS. By tak- 
ing the ratio Igg{G)/ Ipfs{G) at a given reciprocal lattice 
vector G, then this DW factor drops out. Due to the 
presence of this coupling, the (u^) in SS is different than 



ieijqiUj(q^LOn)/q are the longitudinal and transverse 
components of the displacement field respectively. Note 
that Eqn[T2] shows that only the longitudinal component 
couples to the superfluid 9 mode, while the transverse 
component is unaffected by the superfluid mode. This 
is expected, because the superfluid mode is a longitu- 
dinal density mode itself which dose not couple to the 
transverse modes. 

From Ean ll2[ we can identify the longitudinal- 
longitudinal phonon correlation function: 



, ^ ^ KCjg + psq^ 

^""'"'^ (^^2 + p^q^^pn^l -f (A + 2p)q^) + aW^l 

(13) 

The {99) and {ui9) correlation functions can be simi- 
larly written down. By doing the analytical continuation 
iujn LJ + iS^ we can identify the two poles of all the cor- 
relation functions at Wj. = v'i^q^ where the two velocities 
v± is given by: 



(14) 

I 



CO" 




FIG. 3; The elementary low energy excitation inside a ex- 
citonic supersolid (ESS). The coupling between the phonon 
mode ujp = Vpq ( the upper dashed line ) and the superfluid 
mode u!s = VsQ { the lower dashed line ) leads to the two new 
longitudinal modes ll}± = v±q ( solid lines ) in the ESS. The 
transverse mode stays the same as that in the ENS and is not 
shown. 



that in NS, so the DW factor will not drop out in the 
ratio. In this section, we will calculate this ratio and see 
how to take care of this factor when comparing with the 
X-ray scattering data. 

The density order parameter at the reciprocal lat- 
tice vector G is pg{x,T) = e^G-u{x,T)^ then (p^(x, r)) = 



J 



V^. = [k(A + 2/i) + PsPn + 0^ ± sj (k(A -I- 2/z) -I- PsPn + O^Y - ^t^PsPni^y + 2/i)]/2Kp, 
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g-iG,G,{«,«,)^ The Debye-Waller factor: 

/(G) = |(pg(x,r))|2 = ^-G.G,{u.iS,T)u,(x.T)) 

where the phonon-phonon correlation function is: 

(uiUj) = {uiui)qiqj + {utUt){Sij - qiq,) 



where qiqj ^ 



(15) 



(16) 



Then substituting EqnlTB]into Eqn[Tn] leads to: 

- „-iG=[(u?(S,r))ss-(«f(S,r))„s] 



y(G) = Iss{G)/Ins{G) = e-^' 



where the transverse mode drops out, because it stays 
the same in the SS and in the NS. 



Defining {/S.u'^)i{q,iu}n) = {\ui{q,iujn)\'^) ss ~ 

{\ui{q,iuJn)\^)NS, {Au'^)i{q) = J2iujJ'^u^)i{?^^^n) 

and {Au^)i = (uf{x,T))ss - {uf{x,T))NS = 

Ij^^pE^^J^u^W.^^n) = J^iAu'm, it 



(17) is easy to see: 
I 



(fq 1 



R ^ 



9 9 1 

-a q^n 



(27r)2 [3 [{nu:l + p.q^){pnOjl + (A + 2p)q^) + a\^uj^^\p,,u,l + (A + 2p)q^\ 



(18) 



Obviously, (Au^); < 0, namely, the longitudinal vibra- 
tion amplitude in SS is smaller that that in NS. Then 
a(G)(T = 0) = e-^G'(An^), > ^ xhis is expected, be- 
cause the SS state is the ground state at T < Tss, so the 
longitudinal vibration amplitude should be reduced com- 
pared to the corresponding NS with the same parameters 
Pn, \p- 

After evaluating the frequency summation in Eqn ll8[ 
we get: 



(Au^iT) 



(fq 1 coth/3w+q/2 cothf3vpq/2 



(2^)2 p„ 2v+q 



2vpq 



^ Vs-v^^ ^^ cothf3v+q/2 _ coth(3v^q/2 ^^^^ 



V 



+ 



v_ 



2v+q 



2v-q 



At r = 0, the above equation simplifies to: 



{Au')i{T = Q) 



<fq I ^ I 1 
(27r)2 p„ 2v+q 2vpq 



-)( 



1 



1 



__ {v+ + V- - Vp - Vs) A 
{v++v-)vp Airpn 

Z2 1 



A 



(w+ +v- +Vp + Vs){v+ + v-)vp inp, 



< 0(20) 



where A ~ 1/a is the ultra-violet cutoff and we have used 
the fact + i;_ > Vp + Vg. 

At finite T, in EgnfTUl if taking q limit, it is easy 
to see the integrand identically vanishes. So EquHni is 
well defined. 



VII. 



DENSITY-DENSITY CORRELATIONS 



The density-density correlation function in the SS is: 

{p-[x, t)p*~{x' ,t')) = e-kC!iGj{(ui{x,t)-Ui(x' ,t'))(Uj(x,t)-Uj(x\t'))) 

(21) 

where t is the real time. 

For simplicity, we only evaluate the equal-time corre- 
lator {p^{x,t)p*g{x',t)) = {pg{x,T)pg{x' ,t)) where t is 
the imaginary time. It is instructive to compare the den- 
sity order in SS with that in a NS by looking at the ratio 
of the density correlation function in the SS over the NS: 



(f - f ) = {PaP}i)ss/{PaP*^)NS = e-i^'^^^(^-^"> 

(22) 



It is easy to find that 

Ai^p(f-f')=/(0(2 

(23) 

where {Av?')i{q) is defined above EqnfTSland is the inte- 
grand in Eqn IlQI 

At T = 0, the above equation can be simplified to 



^^q■{x-x ))(^y2-)^(^ 



ADp{x-x') 



_ {v+ -Vp- Vs) 1 



{v++V-)Vp 2'Kpn\x—x'\ 



n (24) 



So we conclude that ap{x—x') < 1, namely, the density 
order in SS is weaker than the NS with the corresponding 
parameters Pn, A, p. This is expected because the density 
order in the SS is weakened by the presence of moving 
vacancies. 

It is known that there is only algebraic order instead 
of true solid order at any finite T, so the algebraic decay 
exponent in ESS stays the same as that in ENS. 
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FIG. 4: Specific heat in ENS and ESS. The area below C„/r 
from r = to T = Tss is the entropy in the ESS. The excess 
entropy is due to the lower branch 



VIII. SPECIFIC HEAT IN THE ESS 

It is easy to see that at low T, the specific heat in the 
ENS 



C 



NS 



riNS 



where C^^^ 



r<NS 

3T\2 



Cy 



(25) 



f C(3)A:B(|g^)^ is due to the longitu- 



riNS 
'^tp 



dinal phonon mode where C(3) = 1.202 and 
f C(3)^s(lf7")^ is due to the transverse phonon mode, 
while Cyan is from the vacancy contribution. It is still 
not known how to calculate Cyan yet. The specific heat 
in the SF Cf^ = |C(3)fcB(i^)^ is due to the SF mode 
6. From Eqn[T2l we can find the specific heat in the ESS: 



c 



ss 



TT nV-L 



-C(3)fcs(^)2 + C*^' (26) 



where C*p stands for the contributions from the trans- 
verse phonons which are the same as those in the ENS. 

It was argued in^, the critical regime of finite temper- 
ature NS to SS transition in Fig.l is much narrower than 
the that of SF to the NL transition, so there should be a 
jump in the specific heat dXT — Tss- EgnESl shows that 
at T < Tss, the specific heat still takes ~ T^ behavior 
and is dominated by the uj^ mode in Fig. 2. 

From Eqn l261 it is easy to evaluate the entropy inside 
the ESS: 

Sss{Tss) = ^Cmsi^^n^ + ^) + S'P (27) 



2tt 



where 5*^ is the entropy due to the transverse mode. 

From Eqn l25l and 1261 we find that the excess entropy 
due to the vacancy condensation is: 



AS = 



Tss 



dTCyan/T = |-C(3)A:b(^^^)'(^ + ;^ 



IX. VORTICES IN THE EXCITONIC 
SUPERSOLID 

In section 5, we studied the low energy excitations 
shown in the Fig. 2 by neglecting the topological vortex 
line excitations. Here, we will study how the vortex line 
interaction in ESS differ from that in the ESF. We can 
perform a duality transformation on Egn llll to the vortex 
representation: 



r — 

^y — 



1 



■{efiiyxduax - adaUaSp^r)^ + i2TTaf^j'" (29) 



2K, 



where ^, i^, A stand for space and time, while a,P stand 
for the space components only, Kq = k, Ka = Ps and 
is the 3 component gauge fields and — ^ep^uxdudxO 
is the vortex line current due to the topological phase 
winding in 9. 

Ean l29l has the gauge invariance — > -I- 
where x is any field. It is most convenient to choose 
the Coulomb gauge daaa = to get rid of the longi- 
tudinal component, then the transverse component is 

Q>t — ICapQa 

ap/q. It can be shown that |atp = \aa\'^. 
Then Eqnl^ becomes: 



[Ly 



-^VPn^n 



(A + 2^ + a7K)g2]|u;(g,tJ„)|^ 
2 



+ ^('?V« + ^'M)|at|' + :f- 
2 2p, 

- aq^ / KUi[-q,~uJn)at{q,uJn) 
+ i2nj^a(i + i2TTj^aa 



\ao\ 



(30) 



where the transverse phonon mode ut was omitted, be- 
cause it stays the same as in the NS as shown in EqnfT^ 
It is easy to see that only at has the dynamics, while 
ao is static. This is expected, because although has 
3 non- vanishing components, only the transverse com- 
ponent at has the dynamics which leads to the origi- 
nal gapless superfluid mode uj"^ ~ vlq^ . Eqnl30l shows 
that it is the longitudinal phonon mode ui coupled to 
the transverse gauge mode at- The vortex line den- 
sity is Jq = -^iapdadpO and the vortex current is 
■?a = -^^ai3[do,di3]6. Integrating out the ao, we get 
the vortex density-density interaction: 



y\jUy,T) (31) 



v: 
(28) 

Obviously, AS* > is due to the lower branch v- < 
vip. At a = 0, the above equation reduces to AS" = 
^C(3)fcB(^^^|^)^ which is simply dur to the vacancy 
condensation into ESF. 



Namely, the vortex density- density interaction in SS 
stays as ~ Inr which is the same as that in NS ! The 
vortex energy and the critical temperature Tss in Fig- 2 
^ is solely determined by the superfluid density ps only, in- 
— 2-)iependent of any other parameters in the action Eqn ll2l 
^'p The critical behaviors of the vortices close to the 2d XY 
transition was just the well known KT transition. 

Integrating out the , we get the vortex current- cur- 
rent interaction: 



{2Trf /2j^{-q, -uJn)Daf3iq,uj„)jf!{q,uJn 



(32) 
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where Da^(g,w„) = {S^p - ^||^)A(g,w„) where 

Dt{q,ujn) is the at propagator. Defining ADt{q,u!n) — 
Df {q, u!n) — Df^ {q, a;„) as the difference between the at 
propagator in the ESS and the ESF, then from Eqnl30l 
we can get: 



Namely, the vortex current-current interaction in SS is 
stronger than that in the SF with the same parameters 



X. DETECTION OF THE ESS BY ITS 
PHOTON EMISSION PATTERNS. 

In the previous sections, we worked out the fundamen- 
tal properties of the ESS and various experimental sig- 
natures of the ESS. These experiments can be performed 
easily in the gate voltage generated exciton systems in^iS., 
but it is very hard to perform for the photon pumping 
generated excitons systems inii^i^ due to the short life 
time of these kind of excitons. One important question 
is how to detect the existence of the metastable ESS in 
Fig.l if it indeed exists ? In a recent preprint^, the 
authors found that when the angle resolved power spec- 
trum (ARPS) from the ESF in the Fig.l always takes a 
macroscopic superradiance form which is proportional to 
N'^. The concept of super-radiance was first proposed by 
Dicke in 1954 for N two level atomic atoms confined into 
a small volume V which is much smaller than the wave- 
length of the emitted photoni^i. The macroscopic super- 
radiance from the ESF even holds in the thermodynamic 
limit where ^ oo,y ^ oo, but keep N/V at finite. 
This is due to the macroscopic phase coherence of all the 
N excitons in the whole volume V , so it a completely 
different mechanism than that of Dicke in conventional 
quantum optics. This fact can be used to detect the ESS 
easily in a ARPS experiment for the following reason: the 
power spectrum from the ENS in Fig.l is just a normal 
radiance, so proportional to Njys, however, that from 
the ESF in Fig.l is a macroscopic super-radiance, so pro- 
portional to N^g^^. Because the ENS is a uniform phase 
consisting of a normal solid component and a superfluid 
component, the ARPS at any given in-plane momentum 
k from the ESS is: 

SESsik)<xN^,, + NNs (35) 



For simplicity, we just give the expression for the equal 
time: 



(34) 

I 

which takes a Lorentzian form with the photon energy 
— Vg \Jk'^ + k1 centering around the exciton gap Eg ~ 
1.545eyi6. 

For a simple estimate, the superfluid component is 
only a small ~ 10~^ component of the whole system, so 

N,ac/NNS then Ssssik) (x Nns{10-^Nns + !)• 

Although 10^^ is a very tiny number, but it times a huge 
number N^s ~ 10^° for a 1cm? sample. This leads to 
the fantastic fact that even the ESF component is just 
a very small component of the ESS, it dominates the 
ARPS due to its macroscopic super-radiance. So as dis- 
tance changes in Fig.l, the ARPS along any direction 
should distinguish the three phases ESF, ESS and ENS 
easily and without any ambiguity. 

It was shown ini^, due to the symmetry breaking in 
the ESF, the emitted photon along the normal direction 
is a coherent state. The photon number is proportional 
to the condensation fraction Nq. It was also shown that 
due to the non-vanishing anomalous Green function in 
the ESF state in Fig. la, the emitted photons along all 
tilted directions are always in a two mode squeezed state 
between k and —k even off the resonance. The squeezing 
ratio is proportional to the number of excitons N in the 
EHBL. All the results are robust and independent of any 
microscopic details. The detection of the coherent state 
along the normal direction and the squeezed state along 
all the tilted directions requires phase sensitive homo- 
dyne experiment which is essentially a phase interference 
experiment between the emitted light and a local refer- 
ence oscillator. When increasing the distance as shown 
in the Fig.l, the coherent photons along the normal di- 
rection will be much reduced from the ESF to the ESS by 
10~^ and completely disappear inside the ENS. The two 
modes squeezing ratio will also be much reduced from the 
ESF to ENS by 10~^ and completely disappear inside the 
ENS. However, due to the tiny superfluid component in- 
side the ENS, the phase sensitive homodyne experiment 
to measure the two mode squeezing state inside the ENS 
is much harder to perform than the angle resolved power 
spectrum experiment. 

In principle, the emitted photons from the ESS will 



ADt{x 



X , T 



0) 



AtTk'^P'^ (u+ + f-)(fs + f+)(fs + v^)v+v^Vs \x — x'\ 
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form a photonic band structure reflecting the periodic tri- 
angular lattice structure of the ESS. However as shown 
iit^, the maximum in-plane wavevector of the emitted 
photons around energy Eg ~ 1.5A5eV is at least 2 or- 
der of magnitude smaller than the size of the Brillouin 
Zone of the ESS, so the effect of the underlying lattices 
is negligibly small. 

XI. CONCLUSIONS. 

We studied phases and quantum phase transitions in 
the dilute limit ( 1 ) as distance is increased. When 
the distance is sufficiently small, the system is in the Ex- 
citonic superfluid state. When the distance is sufficiently 
large, it is in the weakly coupled Wigner crystal state. 
We argued that there is a Excitonic supersolid state in- 
tervening between the two phases at some intermediate 
distance range y/r^ < d/as < rg. We derive a quantum 
Ginsburg-Landau action to describe the ESF to the ESS 
transition driven by the collapsing of the roton minimum. 
In general, there should be zero-point quantum fluctua- 
tions generated vacancies whose condensation lead to the 
formation of superfluid mode inside the ESS. We then 
derive the quantum GL action to describe the transition 
from the ENS at the large distance to the ESS at the 
intermediate distance. By RG analysis, we showed that 
the coupling to the quantum fluctuations of the underly- 
ing lattice is marginal, so the ENS to the ESS transition 
is in the same universality class as superfluid to Mott 
insulator transition in a rigid lattice, therefore has exact 
exponents z = 2^ v = 1/2, 77 = with logarithmic cor- 
rections. We showed that the ESS is a triangular lattice. 
We then study the elementary excitations inside the ESS. 
We found that the coupling to quantum lattice phonons 
is very important inside the ESS and leads to two lon- 
gitudinal supersolidon modes uj± = v±q shown in Fig. 3. 
The transverse modes in the ESS stays the same as those 
in the ENS. Detecting the two longitudinal modes, espe- 
cially, the lower branch tj_ mode by neutron scattering 
or acoustic wave attention experiments is a smoking gun 



experiment to prove or disprove the existence of the ESS 
in the intermediate distance regime. Then we calculated 
the experimental signature of the two modes. We found 
that the longitudinal vibration in the ESS is smaller than 
that in the ENS ( with the same corresponding solid pa- 
rameters ), so the Debye- Waller factor at a given recip- 
rocal lattice vector is larger than that in the ENS. The 
density-density correlation function in the ESS is weaker 
than that in the ENS. The specific heat in the SS is still 
given by the sum from the transverse phonons and the 
two longitudinal phonons and still shows ~ behaviors. 
The contribution from longitudinal part is dominated by 
the lower supersolidon branch in Fig. 3. By going the 
to the dual vortex representation, we found the vortex 
density-density interaction in ESS stays the same as that 
in the ESF ( with the same corresponding superfluid pa- 
rameters ), so the vortex energy and the corresponding 
ESS (A) to ENS (A) transition temperature Tkt in Fig. 2 
is solely determined by the superfluid density and inde- 
pendent of any other parameters. The vortex current- 
current interaction in stronger than that in the SF. In 
principle, all these predictions can be tested by experi- 
mental techniques such as X-ray scattering, neutron scat- 
tering, acoustic wave attenuations and heat capacity. 

For the excitons generated by photon pumping, the 
ESS is a meta-stable non-equilibrium state, it will de- 
cay eventually by emitting photons. The angle resolved 
power spectrum ( ARPS ) SESs{k) at any given in-plane 
momentum k is dominated by the macroscopic super- 
radiance from its superfluid component, even the super- 
fluid component is just a very small percentage of the 
the whole system. In the ESF, ESS and ENS phases 
in the Fig.l, the ARPS shows the following sequence: 
Sesf > Sess ^ N^ac > Sens ~ N. This se- 

quence can be measured by a ARPS along any direction 
easily and without any ambiguity. 
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of Knowledge Innovation Program (PKIP) of Chinese 
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